PROCEEDINGS 

°F SCIENCE 



Hadron spectrum of QCD with one quark flavor 



oo 

o 

O 
(N 
> 

o 
O 



> 

o 
d 

00 

o 

> 

X 



Federico Farchioni* Gemot Miinster, Tobias Sudmann, JaVr Wuilloud 

Universitdt Miinster, Institutfiir Theoretische Physik, 
Wilhelm-Klemm-Strasse 9, D-48149 Miinster, Germany 
E-mail: |f archion@uni-muenster . de 



Istvan Montvay 

Deutsches Elektronen-Synchrotron DESY, Notkestr 85, D-22603 Hamburg, Germany 

Enno E. Scholz 

Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 

The latest results of an ongoing project for the lattice simulation of QCD with a single quark flavor 
are presented. The Symanzik tree-level-improved Wilson action is adopted in the gauge sector and 
the (unimproved) Wilson action for the fermion. Results from new simulations with one step of 
Stout-smearing (p = 0.15) in the fermion action are discussed. The one-flavor theory is simulated 
by a polynomial hybrid Monte Carlo algorithm (PHMC) at j3 =4.0 corresponding to a = 0. 1 3 fm, 
on 16^' • 32 and 24^ • 48 lattices; the box-size is L ~ 2.1 fm and L ~ 3.1 fm, respectively. At the 
lightest simulated quark mass the (partially quenched) pion mass is ~ 300 MeV. The masses 
of the lightest bound states are computed, including the flavor singlet scalar and pseudoscalar 
mesons Os and Tjj, the scalar glueball 0++, and the A++ baryon. Relics of SUSY in the mass 
spectrum, expected from a large A^^ orientifold equivalence with the ^=1 supersymmetric Yang- 
Mills theory, are discussed. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Hadmn spectrum of QCD with one quark flavor 



1. Introduction and motivation 



The low-energy dynamics of QCD with a single quark flavor {Nf = 1 QCD) strongly differs 
from the one of the multi-flavor theory since chiral symmetry, and its spontaneous breaking, is 
absent. Nevertheless, several open questions of the physical theory may be better understood in the 
theory with minimal number of fermionic degrees of freedom [|l], For example whether 
setting to zero just one quark mass produces physical effects; in the one-flavor model, due to the 
lack of a chiral symmetry, even the definition of the quark mass is non-trivial [|l]]. Another question 
1^] is whether spontaneous breaking of CP is possible for special combinations of the light quark 
masses (negative quark mass for one flavor). 

Another intriguing aspect of one-flavor QCD, emerging from string theory, is the connec- 
tion with the .yV=\ supersymmetric Yang-Mills theory (SYM), also investigated by our collabora- 
tion [^. The equivalence of the two theories in the bosonic sector can be proved at the planar level 
of an "orientifold" large Nc limit Relics of SUSY are therefore expected in Nf = 1 QCD (with 
Nc = 3) where approximately degenerate scalars with opposite parity should be observed. These 
two particles can be identified in the single-flavor theory with the singlet scalar and pseudoscalar 
mesons which we denote with and r\s- 

As we have argued in [^, some chiral symmetry can be recovered mNf = \ QCD by embed- 
ding the theory in a larger partially quenched (PQ) theory with additional valence quark flavors. 
This allows to build pions and obtain a possible definition of the quark mass by the PCAC relation. 

We present here the latest results of an ongoing computation of the low-lying hadron spectrum 
ofNf = l QCD in the Wilson formulation ||^, obtained from 16^ • 32 and 24^ • 48 lattices at lattice 
spacing a ~ 0.13fm (L = 2.1 and 3.1fm). The lightest simulated pion mass is Mj^ ~ 300 MeV. 



2. Partial quenching 

The definition of a quark mass is not immediate in one-flavor QCD due to the lack of chiral 
symmetry (even in the continuum). For example a PCAC quark mass cannot be defined. More 
generally, the bilinear quark operator qq is not protected against scheme-dependent additive renor- 
malizations and the concept of vanishing quark mass could be devoid of physical meaning |jl|]. 
In our approach we take Ny = 2 valence quarks degenerate with the (single) sea quark. With the 
Ny >2 valence quarks plus the sea quark, all mesons and baryons of QCD can be build, appearing 
in degenerate multiplets due to the exact S\]{Nv + 1) flavor symmetry. A PCAC quark mass can be 
defined in the PQ theory by the non-singlet axial-vector chiral current A"^ 

{d*A+p-) 

Relation ( [2. 1[ ) should be considered here just as a possible definition of the quark mass mNf = \ 
QCD (a possible physical meaning of this definition is not claimed at this stage). 

In [^, validity of the GMOR relation for the pions in the PQ extension of the one-flavor theory 
was confirmed: the masses of the pions can be made to vanish by suitably tuning the bare quark 



mass on the lattice; in this situation the quark mass (2.1) vanishes, too. This scenario is confirmed 
by a study of the chiral Ward identities [^J. 
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Table 1: Summary of the runs: the bar indicates runs with Stout-link in the fermion action (see text). 
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Although pions do not belong to the unitary sector of the theory, their properties and the PCAC 
quark mass can be used for a characterization of the one-flavor theory. In particular the low-energy 
coefficients of the chiral Lagrangian can be extracted by an analysis in PQ chiral perturbation 
theory 



3. Simulation 

The gauge action is discretized by the tree-level improved Symanzik (tlSym) lattice action [^] 
including planar rectangular (1x2) Wilson loops. We apply the (unimproved) Wilson formulation 
in the fermionic sector. Previous results were presented for simulations with the original Wilson 
formulation; with the goal of improving stability of the Monte Carlo evolution at small quark 
masses, we now apply Stout-smeared links in the hopping matrix (one step of isotropic smearing 
with coefficient p^v = P = 0.15). 

The update algorithm is a Polynomial Hybrid Monte Carlo algorithm with a two-step polyno- 
mial approximation (TS-PHMC) |[l0[]. A correction factor C[[/] in the measurement is associated to 
configurations for which the eigenvalues of the (squared Hermitian) fermion matrix lie outside the 
validity range of the polynomial approximation. See [^] for more details on the algorithmic setup. 
The sign of the determinant associated to one (light) Wilson quark can become negative on some 
configurations, even for positive quark masses. Since the sign cannot be taken into account at the 
update level, it must be computed "off-line" and included in the correction factor; the expectation 
value of a generic quantity A is therefore given by 

S[dU]o[U]C[U]A[U] 
^""^ ![dU]a[U]C[U] ' ^'-'^ 

where (7[?7] is the sign of the one-flavor determinant. For the computation of the sign we study the 
(complex) spectrum of the non-Hermitian matrix concentrating on the lowest real eigenvalues: sign 
changes are signaled by negative real eigenvalues. We applied the ARPACK Arnoldi routines on a 
transformed Dirac operator. The (polynomial) transformation was tuned such that the real eigen- 
values are projected outside the ellipsoidal bulk containing the whole eigenvalue spectrum [111 



This allows an efficient computation of the real eigenvalues [12]. 
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Table 2: Results for hadron observables mNj = \ QCD (lattice units). 
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The updating has been performed by applying determinant break-up by a factor of two (that 
is, two half-flavors were considered instead of one). The update sequence consisted of two PHMC 
trajectories followed by an accept-reject step by the second (precise) polynomial approximation. 
The precision of the first polynomial approximation was tuned such that an acceptance of about 
90% was obtained. The same acceptance was required for the Metropolis test at the end of every 
individual PHMC trajectory by tuning the trajectory length (0.4-0.5). This resulted in a high total 
acceptance rate of about 80%. The gauge configurations were stored after every accept-reject step 
by the precise polynomial. Relatively short trajectories were chosen in order to have many con- 
figurations for the glueball mass determination. Optimization of the parameters of PHMC turned 
out to have a substantial impact on the integrated autocorrelation times of the average plaquette. 
Information about the sets of configurations generated up to now can be found in Table |l[ New 
results presented in this Contribution concern runs C, A, 6,624- 



We fix the lattice scale by the Sommer scale parameter ro []13[]. For the conversion into physical 
units we assume the conventional value ro = 0.5 fm. As in ordinary QCD, ro/a grows with the 
hopping parameter k; consistently with a massless scheme the value of ro/c? should be extrapolated 
to the critical value K = where the PCAC quark mass ( ^ ) vanishes. The number of quark 
masses at our disposal are however not sufficient for an extrapolation and we rely on the value at 
the highest K for given j3. We obtain in this way a(3.8) ~ 0.19fm and a(4.0) ~ 0.13fm. This 
corresponds to a roughly constant volume L = 2.1 fm. In order to check finite volume effects, a run 
on a larger 24^ -48 lattice has been started, see Table [l], corresponding to L = 3.1fm. 

We observe that, for a fixed quark mass (in lattice units), the Stout-smearing leaves ro/a es- 
sentially unchanged (as can be seen by comparing runs C and A). On the basis of this observation, 
we will assume approximate equality of the lattice spacing for the j8 = 4.0 runs with and without 
smearing; therefore data from both sets will be included in analyzes at fixed lattice spacing. 



4. Hadron spectrum 



The disconnected diagrams of the r]s and Os correlators (with usual interpolating fields) were 
computed by applying for each configuration 20 stochastic sources with complex Z2 noise and spin 
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Figure 1: The mass of the hghtest physical particles in one-flavor QCD as a function of the squared pion 
mass in units of the Sommer scale. The hadron masses are multiplied by the value of ro/a at the given k. 

dilution. In consideration of autocorrelations, we analyzed every 10th, 16th or 32th configuration 
according to K. The resulting statistics vary between 400 and 600. 

In the case of the baryon, we take the interpolating field A,(x) = EabcWai^Y C'Yiyifb{x)\^lfc{x) . 
The low-lying projected state is expected to be a spin 3/2 parity -positive particle which we denote 
with Ai (the desired spin 3/2 component is extracted by spin-projection). It corresponds to the 
A++(1232) baryon of QCD if the single quark is identified with the w-quark. 

In the PQ sector we measure the pion observables, Mjcfn, nipcAc, in which case correlators 
trivially coincide with the connected parts of the corresponding rj^ correlators; the nucleon mass is 
determined by applying the standard projecting operator 

No smearing was applied for the extraction of the masses. It turns out however that this 
is necessary in the baryon sector where the approach to the asymptotic behavior is slow. The 
optimization of the overlap with the ground state by Jacobi smearing is in plan. 

Results. The results for the hadron observables in lattice units including new runs with Stout- 
smearing are reported in Table ^. The lightest pion mass corresponds to ~ 300MeV. From the 
comparison between runs B and S24 (L = 2.1 fm and L = 3.1fm, respectively) finite volume effects 
can be estimated. These are below the statistical accuracy (~ 10%) for the particles in the unitary 
sector. In the pion sector, they exceed statistical accuracy in the case of the pion mass and of the 
PCAC quark mass, however they are relatively small (~ 4% in both cases). Observe that the sign 
of the finite size scaling agrees with what is observed in standard QCD. 

In Fig. [1] the hadron masses in the unitary sector are reported in physical units as a function of 
the squared pion mass. The a, and the A^ masses are (surprisingly) near to the values observed in 
nature for the corresponding particles or resonances. In contrast, the rjg meson is much lighter than 
the QCD flavor-singlet rj'; this can be understood in terms of the Witten-Veneziano formula. 
The rjs mass, which shows a clear quark mass dependence, can be extrapolated to zero quark 
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Figure 2: The PQ sector of one-flavor QCD: the nucleon mass as a function of the square pion mass in units 
of the Sommer scale. The pion mass is also reported for comparison. 



mass by applying the LO chiral perturbation theory formula [[T^, ^ M^^ = (M^ + a) , 

with and a constants. The result is: roM^^im^ = 0) = 0.84(5) [330(20) MeV]; the error also 
includes a rough estimate of the extrapolation uncertainty, obtained by comparing results including 
or excluding the heaviest quark mass. It is interesting to compare this result for the rj^ mass with 
the one obtained by our collaboration for the corresponding particle in SU(2) SYM, the adjoint rj' 
(a-rj'): r^Ma-r]' = 1.25(5) [499(20) MeV] [Q]. Since the quark mass dependence of the mass 
cannot be seen with our statistical accuracy, we simply apply a fit to a constant; the result is: 
roMf,^ =2.05(9) [810(35)MeV]. 

We can now check our results against the prediction from the orientifold planar equivalence 
Mj^jMo, = {Nc - 2)/Nc X (1 + 5), where the correction 5 = 0{l/Nc,l/N^) is expected to 
be suppressed. We obtain: 



^ =0.410(32)(25), 5 = 0.23(10)(7) 



(4.1) 



(the second error comes from the extrapolation, only data from the L = 2.1 fm volume are included). 
The relatively small value of 5 seems to confirm suppressed 0{\/Nc) corrections, as expected on 



theoretical grounds ||15|]. Small deviations from the leading orientifold prediction were also ob- 
served for the fermion condensate, in numerical simulations [[l6|], and in an analytical computation 



for free staggered fermions [17]. The inclusion in the analysis of smaller quark masses and larger 



volumes could further lower the ratio ( pl| ) and therefore the deviation 5. 

Purely gluonic operators, the glueballs, project onto Spin states, too. We investigate here 
the 0++ state, which is expected to mix with the a,. We neglect for the moment possible mixings 
with the mesonic state and consider diagonal correlators only. Since the computational load is low, 
we analyze in this case each configuration; this allows to obtain a decent signal for the case of the 
Stout-runs, where fluctuations are reduced. We obtain roMo++ = 1.94(25) for run A and 2.43(35) 
for run B; these results are in the ballpark of the Os mass in accordance with strong mixing. 
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The behavior of the nucleon mass as a function of the pion mass squared is reported in Fig. ^; 
for comparison, we also report the pion mass. Also in the case of the nucleon we observe a surpris- 
ing agreement with the expectations from the physical world. 

5. Conclusions and perspectives 

New data from simulations of lighter quark masses in sufficiently large volumes allowed first 
quantitative estimates for the hadron spectrum of Nf = 1 QCD. Results for the Os and T]^ masses 
could be compared with the predictions from the orientifold equivalence; the deviation from the 
leading formula for the ratio of the masses turns out to be relatively small, in accordance with 



observations for the fermion condensate [|1^, |17|]. With the exception of a lighter T]^, no other strik- 
ing deviation from the physical (multi-flavor) picture is observed for the measured quantities. The 
simulation of additional lighter quark masses in view of an analysis in chiral perturbation theory is 
planned for the future. The search for the expected CP- violating phase transition is ongoing. 

The computations were carried out on Blue Gene L/P and JuMP systems at JSC JUlich (Ger- 
many). 

References 

[1] M. Creutz, Annals Phys. 322 (2007) 1518. 
[2] M. Creutz, Phys. Rev. Lett. 92 (2004) 162003. 
[3] M. Creutz, Phys. Rev. Lett. 92 (2004) 201601. 

[4] K. Demmouche, F. Farchioni, A. Ferling, I. Montvay, G. Miinster, E. E. Scholz, J. Wuilloud, 
PoS(LATTICE 2008)061. 

[5] A. Armoni, M. Shifman and G. Veneziano, Nucl. Phys. B 667 (2003) 170; A. Armoni and 

M. Shifman, Nucl. Phys. B 670 (2003) 148; A. Armoni, M. Shifman and G. Veneziano, Phys. Rev. 
Lett. 91 (2003) 191601. 

[6] F. Farchioni, I. Montvay, G. Miinster, E. E. Scholz, T. Sudmann and J. Wuilloud, Eur. Phys. J. C 52, 
305 (2007); PoS(LATTICE 2007)135 [arXiv:07 10.4454 [hep-lat]]. 

[7] C. Bernard, M. Golterman, Y. Shamir and S. R. Sharpe, Phys. Rev D 77 (2008) 1 14504. 

[8] P Weisz, Nucl. Phys. B 212 (1983) 1; P Weisz and R. Wohlert, Nucl. Phys. B 236 (1984) 397. 

[9] C. Morningstar and M. J. Peardon, Phys. Rev D 69 (2004) 054501. 

[10] I. Montvay and E. Scholz, Phys. Lett. B 623 (2005) 73; E. Scholz and I. Montvay, PoS(LAT2006)037. 

[11] Y. Saad, "Numerical Methods for Large Eigenvalue Problems", Manchester 1992. 

[12] H. Neff, Nucl. Phys. Proc. Suppl. 106 (2002) 1055. 

[13] R. Sommer, Nucl. Phys. B 411 (1994) 839. 

[14] C. W. Bernard and M. F L. Golterman, Phys. Rev D 49 (1994) 486. 

[15] A. Armoni and E. Imeroni, Phys. Lett. B 631 (2005) 192. 

[16] T. DeGrand, R. Hoffmann, S. Schaefer and Z. Liu, Phys. Rev. D 74 (2006) 054501. 

[17] A. Armoni, B. Lucini, A. Patella and C. Pica, Phys. Rev D 78 (2008) 045019; A. Armoni, B. Lucini, 
A. Patella and C. Pica, in these proceedings. 



7 



